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We derive a complete Feynman diagram expansion for the elastic form factor and the three- 
body photo and electrodisintegration of the three-body bound state using the covariant spectator 
theory. We show that the equations obtained are fully consistent with bound-state equations and 
the normalization condition previously derived for the covariant three-body bound state, and that 
the results conserve current. 
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I. INTRODUCTION 

The covariant spectator formalism pj has been applied 
successfully to the description of NN scattering and the 
deuteron bound state @, Q, the deuteron form factors 
0, 0] electrodisintegration of the deuteron 0], and the 
three-body system pj, with emphasis on the numerical 
solution of the three-body bound state equations 0, 0] . 
While most of the theory has been developed using di- 
agrammatic methods, all of the two-body theory [lj} 
and the covariant normalization condition for three-body 
spectator wave functions have also been derived al- 
gebraicially. Kvinikhidzc and Blankleider have also stud- 
ied the current for both the two and three-body spectator 
equations [13,1131 using the gauging of equations method. 

In this paper we use the same diagrammatic method 

that is the basis of the results given in Refs. a a 

a a a a a 01 to obtain the complete set of gauge in- 
variant Feynman diagrams needed to evaluate the three- 
body form factors and three-body breakup amplitudes 
for three-body photo or electrodisintegration. Our re- 
sults agree with those recently dervied algebraically by 
Adam and Van Orden 01 1 an d, although they have a dif- 
ferent form, can also be shown to agree with the results of 
Kvinikhidze and Blankleider [l3||. These results, together 
with the relativistic wave functions solutions found in 
Ref. 0, are currently being used to study the three 
body electrodisintegration process 3 He (e, e 'pp)n recently 
measured at the Jefferson Laboratory [lq . 

This paper is divided into five sections and a long Ap- 
pendix. Following this short introduction, we first review 
the definitions and some properties of the covariant three- 
body spectator subamplitudes and spectator equations. 
In Sec. Ill, the basic results for both the elastic and in- 



elastic three-body currents are derived from an analysis 
of the infinite series of Feynman diagrams that consis- 
tently defines both the scattering equation and the cur- 
rents. The derivation is based on the gauging method 
of Kvinikhidze and Blankleider ^| (referred to as KB 
throughout this paper), with a detailed diagrammatic 
proof of gauge invariance and a full discussion of dou- 
ble counting. Then, in Sec. IV, explicit algebraic forms 
of the final results for the currents are presented and in 
Sec. V conclusions, and a detailed comparison with KB 
are given. An algebraic proof that the elastic current con- 
serves the charge of the three-body bound state is given 
in the Appendix. 



II. WORKING WITH SPECTATOR 
THREE-BODY AMPLITUDES 

A. Covariant Faddeev subamplitudes 

In this paper {i,j,k} denote any permutation of par- 
ticles 1, 2, 3, so that j 7^ i, k ^ i and j ^ k, and, for ex- 
ample, i can represent any of the three particles. Then, 
in the nonrelativistic Faddeev theory, the full three-body 
vertex function, which we denote is decomposed into 
three subamplitudes |r l ) which denote that part of the 
vertex in which particle i is the spectator and the other 
two particles (j and k) were the last to interact. The full 
vertex is a sum of the three subvertices 
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(2.1) 



In the covariant spectator theory , the spectator 
is on-shell, and one of the two interacting particles is 
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l : ^Q* x,= ^^Q™ = ^z^Q*" 
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FIG. 1: (Color online) Figure (a) shows the notation for the vertex 
subamplitude TJ,. Lines marked by an X are on mass-shell, and 
the spectator particle is the one that connects to the small dot 
inside the oval. Note that Vi2^\ = CX\> as shown in (b), and 
V3,2^\ = ?^3' as snown m ( c )- 



also on- shell. Hence there are now 6 possible Faddeev 
subvertex functions, denoted by T*-), where i is the (on- 
shell) spectator and j is also on-shell, so that only one 
of the three particles, k in this example, is off-shell. The 
diagrammatic representation of this amplitude is shown 
in Fig. ^a). If the three particles are identical (as we 
assume in this paper) it can be shown [9j, under the 
permutation operator Vij that interchanges particle i and 
j, that 
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(2.2) 



where £ = ± depending on whether or not the particles 
are bosons or fermions [a diagrammatic representation of 
the first two of these equations is given in Fig. [I]. In prac- 
tice, this means that the particles may be freely relabeled 
in Feynman diagrams (as we will discuss below). 

The subvertex functions in the spectator theory each 
satisfy a different constraint, and care must be taken 
to match this constraint to the physics. For example, 
the energy of particle 3 in the subvertex function |r3,) 
is kzo = Mb — E\ — E2, where throughout this paper 

Ei = 



2 + k 4 will always represent a physical en- 
ergy, Mb is the bound state mass, and are the three- 
momenta of the three particles in the three-body cm. 
However, the energy of particle 3 in the subvertex func- 
tions \rl) and |rf) is E 3 . The energy domains of each 
of the subamplitudes is illustrated diagrammatically in 
Fig. |3 Note that, if Mb < 3m, there are three distinct 
domains that do not overlap. 

We now describe briefly how to derive the bound state 
and scattering equations for three identical particles. To 
make the discussion simple and intuitive, we use a dia- 
grammatic approach (similar to that used in the original 
Ref. 0). 




FIG. 2: The three shaded areas are kinematically allowed regions 
of the energies of the three particles, subject to the constraint that 
their total is Mg < 3m, and that two are on-shell. [To read the 
figure, note that each point on the plane defines three unique ener- 
gies. The energy of each particle is the perpendicular distance from 
one of the sides of the large central triangle (positive if "above" the 
side and negative if "below" , where the positive direction is shown 
by the arrows). The geometry of this Dalitz-likc plot insures that 
the sum of the three energies at each point equals the altitude of 
the triangle (chosen to be Mg).\ The on-shell condition fc^o > m 
requires that fc^o lie beyond the line labeled with the number i at 
each end, and the shaded areas labeled by the number i (equal to 
1, 2, or 3) are regions where particle i is off-shell, and j and k have 
energies E > m. Note that these three areas do not overlap. 




FIG. 3: (Color online) Diagrammatic representation of the con- 
struction of the symmetrized three-body scattering subamplitude 
Tjj (represented by the oval) from unsymmetrized subamplitudes 
T", (represented by rounded rectangles). Only the symmetriza- 
tion of the final state is shown. The 6 figures correspond to the 6 
terms in Eq. (TSl as follows: (a) -> 1, (b) -> CV23, (c) -> CV12, 
(d) -» ^23^12, (e) CPi3, and (d) -> ^12^23- 



B. The three-body scattering equation 

For simplicity, this paper addresses only cases in which 
the three-body scattering is a succession of two-body 
scatterings. This means that three-body forces of rel- 
ativistic origin, as defined in Ref. 7], will be neglected. 

We begin with a review of some of the results of 
Ref. [HI. The quantity of primary interest is the sym- 
metrized three-body subamplitude T% , which can be ob- 
tained from the unsymmetrized subamplitudes . Here 
the superscripts i', i denote the on-shell spectators in the 



(a) + (b) 
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FIG. 5: (Color online) The first order scattering can come from 
only two of the six contributions shown in Fig. |3 [diagrams (a) 
and (b)] and hence enters the symmetrized series with a weight 
of jj. In this figure the oval (to the far right) is the symmetrized 
two-body scattering amplitude, and the rounded squares are the 
unsymmetrized two-body scattering amplitudes. 



FIG. 4: (Color online) Terms up to third order in the expansion 
of the unsymmetrized subamplitudes shown in Fig. |!|] 



initial and final state, respectively, and the subcripts j' ,j 
the on-shell interacting particle in the initial and final 
state, respectively. The symmetrized subamplitude can 
be obtained from the unsymmetrized subamplitude by ac- 
tion of the three-particle antisymmetrization projection 
operator 



A 3 = 



^{l + CPl2 + CPl3 + CP23 
+P12P23 + V 23 Vl 2 } 



(2.3) 



normalized to (^3) 2 = -43- The fully symmetrized am- 
plitude is 



1 n' 1/13 1 33' 1/13 ' 



(2.4) 



Expanding out the final state gives, for example, 



T. 
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(2.5) 



corresponding to the six diagrams shown in Fig. EI Each 
of these 6 diagrams generates 6 more terms when the 
initial state is symmetrized, for a total of 6 x 6 = 36 
terms in all. Equation l|2.2[l holds for both the initial 
and final states, and because of this there is really only 
one distinct subamplitude (all others are related to it by 
a phase), which we choose by convention to be T 2 1 2 1 . 

The symmetrization process introduces various weight 
factors into the power series expansion of the amplitude 
T 2 1 2 1 , and hence into the equation for this amplitude. 
These were derived algebraically in Ref. |ll|: here they 
are obtained diagrammatically from a study of Figs.0HH 
Figure 01 shows how each of the unsymmetrized subam- 
plitudes appearing in Fig. is expanded up to third or- 
der (in the two-body scattering). Considering the sym- 
metrization of the final state only, there are a total of 
two first-order terms, four second-order terms, and 8 
third-order terms. The contributions from the first-order 
terms are illustrated in Fig. [SJ the end result is that the 
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FIG. 6: (Color online) Final form of the series for the symmetrized 
subamplitude T^ 1 , with the proper weights for the contributions 
from each order. 



symmetrized first-order scattering diagram enters with a 
weight factor of |. Applying the same argument to the 
other terms gives the series shown in Fig. El which is un- 
changed by the symmetrization of the initial state. This 
series results from the iteration of the scattering equation 
shown diagrammatically in Fig. and hence this is the 
correct scattering equation. This is the result obtained 
previously in Ref. 11]. 

In algebraic form, this scattering equation is 



mil 
J 22 



-Ml 2 - 2tMl 2 G\V 12 T£ . 



(2.6) 



identical to Eq. (3.55) of Ref. |ll|. Here M 22 is the sym- 
metrized amplitude for the two-body scattering of parti- 
cles 2 and 3 [with particle 1 a spectator], and G\ is the 
propagator for particle 3 off-shell. The amplitude M 22 
satisfies the integral equation 



M\ 2 



V 2 \ 
V 2 \ 



V 22 G\ M 22 

^22 ^2 ^22 ' 



(2.7) 
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FIG. 7: (Color online) Diagrammatic representation of the equa- 
tion for the spectator three-body scattering subamplitude T22 ■ 
This equation is equivalent to the series shown in Fig. |g] 



FIG. 8: (Color online) Diagrammatic representation of the equa- 
tion for the symmetrized two-body scattering subamplitude M^- 
The open rectangle represents the irreducable kernel V^. 




(C) 

FIG. 9: (Color online) Diagrammatic representation of the defini- 
tions of three-body subamplitudes with two particles off-shell, (a) 
The scattering amplitude with both particles 2 and 3 off-shell, (b) 
the bound state vertex function with particles 2 and 3 off shell, and 
(c) the bound state vertex function with the spectator (particle 1) 
off shell. 



where is the kernel, or driving terms, of the NN inter- 
action. This equation is illustrated in Fig. |H] In diagrams 
drawn in this paper, the minus sign in the second term 
of (|2.6I) and 1|2.7[) will be associated with the propagator 
G?2, so in the figures the propagator is — G\. However, 
the factors | and 2£ will be shown explicitly. 

The three-body bound state produces a pole in the s — 
P 2 channel (with P M the total momentum four- vector), 

^ = -jjf^l+^ (2-8) 

where the remainder term, TZ, is regular at the pole, and 
the spin structure of the propagating bound state is in- 
cluded in the vertex function jr^). Substituting (|2.8|) 
into |2j3), approaching the pole, and equating residues 
gives the three-body bound state equation 

\Tl) = -2CMl 2 GlV 12 \T\) . (2.9) 

Note that the inhomogenous term in Eq. (|2.6|l has no 
three-body bound state pole, and hence does not con- 
tribute to the bound state equation. The bound state 
equation is therefore diagrammatically identical to Fig. [7| 
without the inhomogenous term. 

We conclude this section with a discussion of how to 
define the spectator amplitudes when one of the spec- 
tators is off shell. The definitions needed in the subse- 
quent discussion are shown in Fig. [!|] Here the princi- 
ple is to expose the two-body interaction which connects 
to the off-shell particle, because the two-body amplitude 
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(A) (a) (b) (c) 




(d) (d) 



FIG. 10: (Color online) The upper panel shows the Bethe-Salpeter 
vertex function as a sum of three subamplitudes for the three dif- 
ferent choices of the last interacting pair. Lower panel shows that 
naive use of this vertex function in the form factor gives overlap 
terms (a) X (b) and (b) X (a) leading (after use of the equation) to 
two terms of type (d), where only one should be present. 



can be extended off-shell by pulling out the last two- 
body interaction (which is always defined with both par- 
ticles off shell). Note that, when the spectator is off-shell 
[Fig.|5{c)], the equation must be iterated twice to get the 
desired result. 

We now turn to the main subject of this paper, the 
diagrammatic derivation of the three-body current oper- 
ator. 



III. SPECTATOR THREE-BODY CURRENTS 
A. The problem of double counting 

To expose one of the central issues in the construction 
of three-body currents, we begin by looking at what ap- 
pears to be the lowest order result in the Bethe-Salpeter 
formalism, and show that this expected result leads to 
over counting. 

The full BS three-body vertex function is the sum of 
three subamplitudes, as shown in Fig. HUT A) . Guided by 
nonrelativistic theory, we might expect the impulse ap- 
proximation to the current to be related to the square 
of the wave function, as illustrated in Fig. HOf B). How- 
ever, if this proposed current is expanded using the wave 
equations, it leads to two terms of type Fig. HOf d). while 
direct examination of the ladder sum (for example) shows 
that there should be only one such term. Unless an in- 
teraction term of type (d) is explicitly subtracted from 
the "impulse" approximation, it will be double counted. 
The same problem does not arise in nonrelativistic the- 
ory because there the diagrams represent a sequence of 
operators which, in general, do not commute. The iter- 
ation of (a) x (b) gives a different contribution from that 
of (b)x (a), and both must be present. [The treatment of 
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(a) (tH) (b 2 ) 

FIG. 11: (Color online) Diagrammatic representation of the gaug- 
ing of the propagator G\ as given in Eq. 13.21 , Diagrams (a), (bi), 
and (D2) are the three terms on the r.h.s. of the equation, respec- 
tively. 

this problem in the context of the Bethe-Salpeter theory 
is discussed in Ref. [l7j-] 

It turns out that the spectator theory, like nonrela- 
tivistic theory, also does not suffer from double counting. 
Furthermore, the topology of the terms shown in Fig. 1101 
can be used to simplify the spectator formalism. A de- 
tailed demonstration of this is best left until after the 
general results have been obtained below. 

B. Step I: Coupling photons to internal lines and 
vertices 

The current is constructed by coupling the photon to 
all propagators and momentum dependent couplings in 
every diagram in the infinite series. This will give, ac- 
cording to a general argument developed by Feynman, 



In each term the particle with momentum p 2 is on shell, 
p 2 — m 2 , and pf = pi ± q, and 

G( P ) = (m- f)- 1 . (3.3) 

Each term is a direct product of Dirac operators on the 
space of particle 2 and 3 (with the space on which the 
operators act implied by the momentum labels, so, for 
example, p 2 operates on the space of particle 2). This 
equation is illustrated diagrammatically in Fig. 1111 

Using Fig.lTTl Eq. [pHl is represented in Fig.[T21 This 
is an equation for (T^ 1 )^, and can be solved by iteration, 
using the series representation for the three-body scat- 
tering amplitude given in Fig. |S| The solution is shown 
in Fig. 1131 To obtain this solution diagrammatically is 
straightforward, if not familiar to many. Fig. 1141 demon- 
strates diagrammatically how Fig. Il^jf c) is obtained by 



a conserved current. After the current has been con- 
structed in this way, a diagrammatic proof of current 
conservation will be given. 

The construction of the current will be carried out in 
three steps. First the coupling to internal lines and ver- 
tices will be constructed, Then the wave equations will 
be used to rearrange the result into a more usable form. 
Finally, the extension of the result to inelastic processes 
requires coupling to the final state nucleons, and the cor- 
rect way to do this will be developed last. 

The coupling to internal lines and vertices is very 
nicely obtained using the gauging of equations method 
of Kvinikhidze and Blankleider The first step in 

this construction is to note that the coupling of the pho- 
ton satisfies the distributive rule of differential calculus. 
Starting from the scattering equation H2.6(l . the photon 
coupling must satisfy the equation 

(ThY = \ (Ml 2 y-2C(Ml 2 y G \V 12 T£ 

-<ltMl 2 G\V^{T},l)\ (3.1) 

where X^ denotes the coupling of the photon to all inter- 
nal lines and vertices in the series of Feynman diagrams 
that make up X. The gauging of the spectator propaga- 
tor generates three terms with the operators connecting 
each of these terms having different arguments. Denoting 
these operators by A and B, the gauged G\ is 



(3.2) 

I 

iterating the inhomogenous term in Fig. I12f c). 

In subsequent applications the discussion will be lim- 
ited to those cases in which the initial state is bound. 
These diagrams are extracted from the general result 
shown in Fig. ED by approaching the bound state pole 
in the initial scattering and retaining the residue. Dia- 
grams (a) and (ai) do not have such a pole, and therefore 
do not contribute. The result for the bound state internal 
current is shown in Fig. 1151 

The parts of this figure involving rescattering in the fi- 
nal state are identical to Fig. [4] of KB [13|. To show this, 
first compare our two-body scattering Eq. I2.7fl with the 
KB two-body scattering equation (Eq. (26) of Ref. [l3|l- 
Note that G\ = —5 2 d 3 , so that the equations are identical 
if 2M 22 = t\ and 2V 22 = v±, corresponding to a different 
normalization of the two-body amplitudes. Next, note 



A(G\YB = 4(p2,p£)G(p+)#(p+ p 3 )G(p 3 )(m+ji 2 )B{p 2 ,p 3 ) 
+A(p 2 ,p 3 ) {m+ i>i) j%(p 2 ,p 2 ) G(p 2 ) G{p 3 ) B(p 2 ,p 3 ) 
+A(p+, P3 ) G(p+)j2(p+,p 2 ) G(p 3 )(m+ i> 2 ) B(p 2 , P3 ) . 
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(a) (b) (c) 




(di) (d 2 ) (e) 



FIG. 12: (Color online) Diagrammatic representation of Eq. 13,11 , 
with (a) the inhomogeneous term (-^22) ' 3 ) ^he t erm 
(m^A' 1 G\Vi2T^ , (c), (di) and (d 2 ) the three terms result- 
ing from (G3;) M as illustrated in Fig. 1111 and (e) the term 




W (d 2 ) 



FIG. 13: (Color online) Diagrammatic representation of the solu- 
tion to Eq. 13.11 for the current arrising from the coupling of the 
photon to all internal lines and vertices. Diagrams (ai) and (t>i) 
include effects of rescattering to all orders arrising from diagrams 
(a) and (b). Similarly, diagrams (c), (di), and (d 2 ) include higher 
order effects arrising from diagrams (c), (di), and (d 2 ) of Fig. 1121 

that our three-body scattering Eq. i|2.6[l is identical to 
the similar Eq. (15) of KB if we choose £ = — 1 and set 
6T22 = Xt%, which corresponds to a different normaliza- 
tion of three-body scattering. With these replacements, 
our diagrams 1151 (bi). (c), and (d) are identical to Fig. 
[4] in KB [recalling that second figure in the KB Fig. [4] 
is the same as the sum of our (di) and (da)]. Since the 
KB result shown in Fig. [4] applies only to transitions to 
connected final states, diagrams like our Fig. ^| (b) will 
be discussed after we have finished our discussion of the 
breakup process. 

Note that Fig. ll5f d9) includes a contribution in which 
the spectator (particle 1) is off-shell. While recognizing 
that further reductions are possible, KB elect to leave 




FIG. 14: (Color online) Figure showing how the the inhomoge- 
neous term, Fig. I12f c'). is iterated by final state interactions to all 
orders, leading to Fig. I13f c1. 





to W «y 



FIG. 15: (Color online) Diagrammatic representation of the so- 
lution, Fig. 1131 applied to processes with an incoming three-body 
bound state. 

their answer in this form. In Step II of our derivation, 
we use the wave equations (|2.Gll and 1)2. 9|) to replace this 
amplitude by an equivalent one in which the spectator 
is on-shell. This replacement does not change the total 
result, but is still very useful for numerical applications, 
and leads to a nice demonstration of how the spectator 
equations avoid the double counting problem in a natural 
way. 



C. Step II: Removal of off-shell spectator 
contributions 

The diagram (d.2) of Fig. ED can be transformed as 
shown in Fig. ^j] The first step is to replace the initial 
bound state amplitude (which has its spectator off-shell) 
using the definition shown in Fig. H2c). Then we recog- 
nize that, if one of the two NN scattering amplitudes 
introduced by this substitution is identified with the fi- 
nal state, and one with the initial state (as outlined by 




FIG. 16: (Color online) Diagrammatic representation of the trans- 
formation of diagram (012) in Fig, 1151 Ths first step uses Fig. J§tc) 
to replace the initial three body amplitude [denoted by (a) in the 
figure]. Then the two parts labeled (b) and (c) are isolated, and 
Fig. EJa) is used to replace (b) and Fig. Elb) to replace (c). The 
final result is the two terms shown in the lower panel. Note that 
the diagram (b2)(c) cancels an identical diagram that comes from 
the first iteration in the final state scattering amplitude in (bi)(c) 
[c.f. diagram (ai) of Fig. 151. 





(b n ) 
(d 2 ) 



(b,) 








FIG. 18: (Color online) Diagrammatic representation of contri- 
butions from internal photon couplings to the three-body breakup 
process. Diagram (a) is constructed from the core diagrams shown 
in Fig. ll9l ('simDlv add the two external half ellipses to each diagram 
in Fig. H9l to complete the diagram). Diagram (b) is the same dia- 
gram (b2)(c) that appeared in Fig. I16l when the off-shell spectator 
was replaced by an on-shell spectator. It cancels the first iteration 
of the final state interaction in Fig. ll9f e9) so that the sum of all of 
these contributions contains no couplings to external nucleons. 



(e) 




FIG. 17: (Color online) Diagrammatic representation of the 
gauged (Mi*)*" (where, in our application, the on-shell particle is 
#2 and the off-shell particle #3). Figures (bi) and (bi) arise from 
the gauging of the propagator (G* ) M , and (c)-(e) include contri- 
butions from the two-body interaction current, represented by the 
photon coupling inside a shaded rectangle. There are no three-body 
forces, and hence no three-body interaction currents. 



the boxes in the figure), it is possible to use the equa- 
tions illustrated in Fig. EJa) and (b) to further simplify 
the diagram. The result is the two diagrams shown in 
the bottom panel. The conclusion is that diagram (d2) 
can be replaced by these two diagrams, neither of which 
has the spectator off-shell. 

It is convenient to further simplify Fig.^|by replacing 
the interaction of the photon internal to the NN ampli- 
tude with the gauged result for the TV TV amplitude, 



FIG. 19: (Color online) This figure shows "core" diagrams com- 
mon to all interactions. The half circles in the initial and final state 
can be either the bound state, or part of a three-body scattering 
amplitude. Substituting into Fig. 1181 gives the complete internal 
current. (A) the core equals the sum of (b) and (d) photon cou- 
pling to off shell particle, either with or without exchange of the 2 
on-shell particles, (c) and (e) the coupling of the photon to the on- 
shell interacting particle (arranged so that the spectator is always 
on-shell), and (f) two-body interaction current diagrams. 





{Ml 2 Y = {V^ + {V^G\Ml 2 



M 22 G\ (^22)1* 



FIG. 20: (Color online) The three-body bound state form factor 
is constructed only from the core diagrams of Fig. 1191 



+M 2 X 2 G\ (V^ G\ M\ 2 + Ml 2 (G\Y Ml 2 , (3.4) 

illustrated in Fig. El This replacement exposes the two- 
body interaction current, (^22)1 which includes the pho- 
ton coupling to all exchanged mesons and meson-nucleon 
vertices (when present). Using this expansion, and the 
wave equations, gives the diagrams shown in Fig. 1 181 with 
the "core" diagrams defined in Fig. ^jj The core dia- 
grams also define the three-body form factor. To extract 
this form factor from the result shown in Fig. ^0 go to 
the bound-state pole in the final scattering amplitude, 
and note that the diagrams (b), (c) and (d) will not con- 
tribute. The result is summarized in Fig. 1201 

It remains now to (i) demonstrate that the form fac- 
tor conserves current, (ii) find the additional diagrams 
(describing couplings to external nucleons) that must be 
added to Fig.^|to compete the description of the three- 
body breakup process, and (iii) show that the three-body 
breakup also conserves current. 

It turns out that the best way to proceed with these 
remaining tasks is to first derive the WT identity for the 
core diagrams, Fig. 1191 
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(a) (b) (c) 



(a^ (a 2 ) (b) (c) 




FIG. 21: (Color online) Diagrammatic representation of Eq. 13.61 . 
The small open triangle on each amplitude marks where the photon 
was inserted, so the momentum of the amplitude at the location 
of the triangle does not match the momentum of the propagator 
connected to it. The larger triangle denotes the sum of two con- 
tributions: (d) is the sum of (ai) and (b); (e) the sum of (a2) and 
(c). 



D. Step III: The WT identity for the core diagrams 

Our study of current conservation is based on a gener- 
alization of the arguments of Gross and Riska ^3 • There 
it was shown that the current will be conserved if (i) it is 
constructed from elementary nucleon and interaction cur- 
rents that satisfy the appropriate WT identities, and (ii) 
contributions from all possible couplings of the photon 
to nucleons and interactions are included in a consistent 
manner. The core diagrams derived in the previous sub- 
section provide a consistent scheme for coupling photons 
to all two-body interactions and internal nucleons, and 
hence provide the solution to condition (ii) when there 
are no external free nucleons (true for the form factor, 
but not for the three-body breakup). It remains now to 
show explicitly how the WT identities for the elementary 
nucleon and interaction currents insure that this is so. 

The nucleon current is constructed to satisfy the WT 
identity 

q^(p',p)= [G- 1 (p)-G- 1 (p')] , (3.5) 

where the particle charge is excluded from the definition 
of the current, and conservation of four- momentum at 
every vertex implies that q = p' — p. From this it follows 
that the divergence of Eq. I|3.2|l is 

q»A(G 1 2 ) tl B = 

A( P2 ,pt) (m+fc)G(pt)B( P 2,P3) 
-A(p 2 ,pi) (m+ fc)G(p 3 ) B(p 2 ,p 3 ) 
+A(p 2 ,Ps) (m+ jfe)G(ps) B(p 2: p 3 ) 
~A(p+, P3 ) G( P3 )(m+ i> 2 ) B{p 2 ,p 3 ) . (3.6) 

This relation is illustrated diagrammatically in Fig. 1211 
Note that Fig. I21f d1 is a shorthand notation for the 1st 



■ - 3p 

FIG. 22: (Color online) Diagrammatic representation of Eq. 13.71 . 
The meaning of the open triangles is explained in Fig. 1211 

and 3rd terms on the r.h.s. of Eq. (|3.6|l and Fig. ^lf e) for 

the 2nd and 4th terms. Similarily, the WT identity sat- 
isfied by the interaction current (for further discussion, 
see Eq. (3.3) of [3) is 

% i v 22 ) i= [V22 ip'2 ,Pi;P2,p 3 )~ V 2 \ {p' 2 , p' s ; p 2 , pf)] 
+ [V 2 \(P2~^P^P2,P3) - Vwip^P^PiiPa)] > ( 3 - 7 ) 

with pf defined above. In this equation, the photon mo- 
mentum is inserted wherever there is a pf , as in Eq. (13.011 . 
Using the notation of Fig. this equation is illustrated 
in Fig. 03 

The derivation of the WT identity for the core dia- 
grams depends on the observation that the three-body 
bound state and scattering equations can be used to ex- 
press the amplitudes in terms of the two-body kernel V 22 , 
as illustrated in Fig. [321 Using this relation, and the WT 
identities (|3.6|l and (|3.7I) . the WT identity for the core 
diagrams can be found. The steps are outlined in Fig.l2"4l 
Cancellations occur when the identitities of Fig. [531 arc 
used to reexpress the diagrams !24f b')- l2~4T e'). In detail, 
Bla) cancels IMe) and IMil. 1Mb) cancels IMf ) and 
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(g) 



(e) 
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FIG. 23: (Color online) Diagrammatic derivation of the replace- 
ment of the scattering amplitude by expressions involving the inter- 
action kernel V 2 2- The top line shows the wave equation with the 
amplitude (A) equal to the inhomogeneous part (a) and the rescat- 
tcring part (b). In the second line, Mli m (b) i s replaced by its 
scattering equation, Fig- |H1 an< I m the third line one of these terms, 
(c), is replaced by the three-body scattering equation. Because of 
the cancellation of diagrams (a) and (f), the final result (last line) 
is the sum of only three diagrams, (g), (e) and (d). The diagrams 
(a), (f), and (g) do not contribute to bound states, leaving only 
diagrams (e) and (d). 




FIG. 25: (Color online) Derivation of the WT identity for the inter- 
nal photon couplings to the three-body breakup process [Fig. 118^ . 
Using Fig. 1241 the divergence of core diagrams leaves diagrams (m) 
and (n). The divergence of llSf bt gives diagram (a), the divergence 
of ll8l c) gives (b) and (c), and the divergence of diagram !18l d) gives 
(d) and (e). Diagram (m) is cancelled by (b) and diagram (n) is 
cancelled by (d), leaving the three diagrams (a), (c), and (e). 



IM j) leaving IMm).Elc) cancels E^g) andEHk), and 
IMd) cancels IMh) and 12117) leaving Efln). 

This result implies that the bound-state form factor 
conserves current, because in this case neither of the di- 
agrams |5Hm) or l^tT n) are present, and the WT identity 
gives zero. Discussion of the three-body breakup pro- 
cess requires additional diagrams, which will be discussed 
now. 



(a) 



-3 



(e) 



(b) 




(i) 



+3 



+6? 




(f) 




+6? 



-12 




(c) 




(g) 




-65 



+ 12 



+6? 




(d) 




(h) 




(i) 



(k) 



(I) 




(m) 



(n) 



FIG. 24: (Color online) Derivation of the WT identity for the 
core diagrams when connected to an initial bound state and a fi- 
nal scattering state. Using Fig. 1211 the divergence of core diagrams 
E0h)-EH e ) gives diagrams (a) through (d) above, and using Fig. 1221 
the divergence of the interaction current diagrams I19f fl gives dia- 
grams (e) through (I). These diagrams reduce to (m) and (n), as 
discussed in the text. 



E. Step IV: Photon coupling to external nucleons 

Using the results of Fig. [2U we obtain the WT iden- 
tity for the internal photon couplings to the three-body 
breakup process, shown in Fig. 1181 The result is shown 
in Fig. H] 

The coupling to external nucleons will produce the 
terms needed to cancel the diagrams in Fig. I25f at . HR c) . 
andESIe). SinceUJa) results only from Fig.nHrbhUSTb) 
will be removed from the final three-body breakup cur- 
rent. [Removal of this diagram means that the core con- 
tributions will contain a diagram like Fig. llSf b') that can 
be interpreted as an interaction with a free final state par- 
ticle, and care must be taken not to overlook this term. 
We will discuss this further in the conclusions.] The final 
result for three-body breakup, including photon coupling 
to external nucleons, is shown in Fig. 1261 

The three-body breakup diagrams (ai) and (a2) will 
be referred to as the relativistic impulse approximation 
(RIA), diagrams (c) and (d) as interaction currents (de- 
noted by I, but not to be confused with isobar currents 
which, if present, are included in the interaction cur- 
rents), and the core contribution (b), which includes final 
state interactions driven by both the RIA and the I (and 
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(c) 



(d) 



FIG. 26: (Color online) Diagrams that describe the three-body 
breakup process. Diagrams (a) are the RIA, (b) the FSI, and (c) 
and (d) are interaction currents (I). 




FIG. 27: (Color online) Proof that the breakup diagrams conserve 
current. Diagram (a) results from application of the WT identity 
13.61 on the RIA terms (with the final state on shell), (b) and (c) 
are from Fig. 1241 and (d) and (e) result from the the application 
of Fig. 1221 on the I diagrams I26f c~) and 1261 d^, respectively. After 
the cancellations shown, the application of Fig. !23l to diagram (a) 
cancels (d2) and (e2), giving zero. 



denoted FSI) . A proof that this set of diagrams conserves 
current is given in Fig. 1271 

We have completed our derivation of the current, and 
return now to the discussion of double counting first in- 
troduced in the last section. 



F. Conclusion: Removal of double counting in the 
Spectator theory 

We now demonstrate that the problem of double count- 
ing referred to in subsection A above is solved by the 
current operator given in Fig. [213 

Guided by Fig. 1101 we examine the "exchange" terms 
I19f e) that contribute to the core process. These dia- 
grams are reproduced in Figs. ^f a - ) and (b) for the case 
when the initial state is bound and the final state is three- 




FIG. 28: (Color online) Diagrams (ei) and (e2) of Fig. ll9l fredrawn 
here as (a) and (b) and without their overall factor of 6£] are equal 
to (c) + (d) + (e). The last line redraws (d) and (e) so that they 
can be easily compared with the diagrams shown in Fig. 1 1 ( ) I 




FIG. 29: (Color online) Fig. 1281 c) (multiplied by the missing factor 
of 6£) is redrawn showing how it could be interpreted in a different 
way. 



body breakup. Using the bound state wave equation with 
Fig. EEJa) and the scattering equation with Fig. I28f b) 
gives the three contributions (c), (d), and (e). Figures 
(d) and (e) are untangled in the last line of the figure, so 
that they may be more easily compared with diagrams 
HUf d) of the Bcthe Salpeter theory. Note that in the 
spectator theory, both diagrams (d) and (e) must occur, 
because they describe different processes, with the spec- 
tator on shell either "before" or "after" the interaction 
(note that here "before" or "after" refer to a topolog- 
ical ordering and not a time ordering). Our insistence 
that the spectator always be on shell has eliminated the 
double counting problem. 

Finally, look at diagram I28f ch which arises from final 
state interactions. In Fig. 1291 this diagram is rearranged 
to look like an RIA contribution with the spectator off 
shell. This would lead to double counting if we had in- 
cluded such processes in the RIA, but these contributions 
are explicitly excluded from the RIA contributions (ai) 
and (a 2 ) shown in Fig. 12^1 

We now record, for future use, the algebraic espressions 
corresponding to our major results. 
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FIG. 30: (Color online) Diagrammatic representation of the equa- 
tion for the symmetrized two-body scattering subamplitude Mg 2 
with both final particles off shell. 




FIG. 31: (Color online) Diagrammatic representation of the equa- 
tion for the three-body vertex function with both interacting par- 
ticles off-shell. 



IV. ALGEBRAIC EXPRESSIONS FOR THE 
WAVE FUNCTIONS AND CURRENTS 

In this section we record the algebraic form of the 
three-body vertex functions, wave functions, form fac- 
tors, and three-body breakup. The Appendix shows in 
detail how the covariant three-body normalization con- 
dition leads to the conservation of charge. 

A. The three-body vertex function and wave 
function 

The three-body vertex function defined in Ref. Q 
[Eq. (3.14)] is denoted 

F\ 1 \ 2a (k 1 ,k 2 ,k 3 ) = (k 1 \ 1 (k 2 X 2 fc 3 a)|r^) , (4.1) 

where (ki\i(k 2 X 2 fc3a)|r|^ is the three-body vertex func- 
tion describing the coupling of the 3 He nucleus to the 
three-nucleon system, with the first pair of arguments 
[fci, Ai] the four-momentum and helicity of the on-shell 
spectator, the second pair [fc2,A2] the four-momentum 
and helicity of the on-shell particle in the interacting pair, 
and the last pair [k 3 , a] the four- momenta and Dirac in- 
dex of the off-shell particle in the interacting pair. Differ- 
ent notations will sometimes be used for these momenta; 
it is their location in the argument list that identifies 
them as spector, on-shell interacting particle, or off-shell 
particle. In this notation the momenta and spin indicics 
of all three particles are defined in the same system, with 

P = fci + k 2 + k 3 , (4.2) 



the total four momentum of the bound state. In the rest 
frame of three-body bound state, the energy of the off- 
shell particle is 

k% = M B -E 1 -E 2 <m. (4.3) 
In this representation, the symmetry of the amplitude is 
simply 

^x 1 x 2 a(ki,k 2 ,k 3 ) = C Vi2^x 1 x 2 a{ki, k 2 , k 3 ) 

= C^x 2 x ia (k 2 ,k u k 3 ). (4.4) 

When solving the three-body bound state equations 
(|2.9() it is convenient to work in the rest frame of the 
bound state, but to boost the interacting pair to its own 
rest frame, where the partial wave decomposition of the 
two-body amplitude that drives the equation, M 22 , is 
defined. The numerical solutions reported in Ref. @ 
were carried out in this mixed frame, i.e., with k\ and 
Ai defined in the three-body rest system and the remain- 
ing variables defined in the rest system of the interacting 
23 pair. The connection between these two representa- 
tions will be discussed in a subsequent paper, where we 
will calculate the electrodisintegration of 3 He [l5j . In this 
section we will use the representation (|4.1I) . 

The wave function, when needed, is defined by 

*A!A 2 a(fci, ^2,^3) = G a ^(k 3 )T Xl x 2al (ki 1 k 2l k 3 ) , (4.5) 
where the nucleon propagator, G, was given in Eq. (|3.3|l . 



B. The three-body form factors 

The diagrams needed to calculate the form factors were 
displayed in Figs. EH and [201 Following work on the 
deuteron form factors [5j, the diagrams Figs. EJb)-(e) 
are referred to as the complete impulse approximation 
(CIA). Diagrams Figs. lliJf f) are the interaction currents, 
denoted by I. 

Some of the CIA diagrams require knowledge of the 
vertex function with the two interacting nucleons off- 
shell. This vertex function was defined in Fig. Elb). A 
more convenient expression for this vertex function can 
be found using the equation for the two-body scattering 
amplitude with both particles in the final state off-shell 
[this generalization of Fig.[S]is shown in Fig.[3U|. Substi- 
tuting Fig.|2n]into Fig-EIb), and using Fig.^b) a second 
time, gives the result shown in Fig.|^for the bound-state 
amplitude with particles 2 and 3 off-shell. [This is a gen- 
eralization of a result previously shown in Fig. 1231 ] The 
algebraic form of the equation shown in Fig. 1311 is 



Fx^aiki, k 2 , k 3 ) = - 



m d 3 k' 2 \ ~\ 



V, 



(fc 2 ,fc s 



;fc^ 3 )[l + 2CPi2] Vx^cfaXX): 



(4.6) 
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where V is the same two-body kernel used in Eq. (|2.7() . and summation over repeated Dirac indicies is implied. The 
phase of each spectator amplitude is computed from £(—£)", where n is the total number of off-shell propagators plus 
interactions (either vertices or kernels) in the integrand of the amplitude. Equation Ij4.6|l has one off-shell propagator 
(contained in 'J), and two interactions (V and T), for a phase of i(— i) 3 = —1. For convenience we have adopted the 
notation 

V f 3 a .f3' a '{k2 1 k 3 ;k' 2l k' 3 )u, }l (k' 2 ,X' 2 ) = V /3at x' 2a '(k 2l k 3 ;k 2 ,k 3 ) 

Ufs(k 2 ,X2)Vp aM a'(k2,k3;k' 2 ,k 3 ) = Vx 2a ,x' 2 a' (k 2 , h; k' 2 , k' 3 ) , (4.7) 

so care must be taken to distinguish Dirac indicies from hclicity indices. Whenever a Dirac index is replaced by a 
helicity index, a contraction with an on-shell, positive energy spinor, such is shown in Eq. I|4.7II . is implied. The 
on-shell Dirac spinors are normalized to fin = 1. 

The algebraic result for the six diagrams that make up the CIA can now be written 

^ia = 3e|| ^^ity E {^w(fci, fa, fc 3 + ) [1 + 2 CP12] & a (kl,k 3 ) * AlA2a (fcx, fa, fa) 

Ai A2 

+Tx 1 p' a (fa,k 2 h , k 3 ) Gp> p{k%) ,fa) [1 + 2CPi 2 ] u 7 (fa, A 2 ) ^x 1 x 2 a{fa,fa,fa) 

+*A 1 A 2 c e (fci,fe,fc 3 f Wfe,A 2 )[l + 2C^ (4.8) 

where the doubly off-shell vertex functions are evaluated using 14. j a > a (k',k) is the single nucleon current for off- 
shell nucleons with incoming (outgoing) four-momenta k^ (k'^), — k% ± q, and in every term k\ = fcf = m 2 and 
ki + k 2 + k 3 = P, where P is the four-momenta of the incoming deuteron. Each diagram has two off-shell propagators 
and three interactions, for a phase of i(— i) 5 = 1. The I diagrams are 

TtI „ fff m 3 d 3 kid 3 k2d 3 k! ? v-^ - ., ,, ,, w , „ ^ , 

Jl = 3e JJJ Ek&EttLy E*^h,k' 2 ,k>)[l + 2tV 12 ] 

Xi A2 A 2 

x ^, Q , ;A2 , Q (fc 2 ^3;fe,fc 3 )[l + 2C7'i2]*A 1 A 2Q (A ; i,A :2 ,fc 3 ), (4.9) 
I 



where fa + fe 2 + fa = P, fa + k' 2 + k 3 = P' = P + q, 
kj = k' 2 2 = k 2 =m 2 , and V^^.^ Jk!,,^; fa, fa) is the 
symmetrized two-body interaction current for nucleons 
with incoming four- momenta fa , k 3 and outgoing four- 
momenta k' 2 ,k 3 . In Vi the first four-momentum listed 
in each pair describes an on-shell nucleon with incoming 
helicity A 2 and outgoing helicity A2 . The second nucleon 
is off-shell, with incoming Dirac index a and outgoing 
Dirac index a' . 

The equations l|4.8|) and l|4.9l) . with the two particle 
off-shell vertex function defined by Eq. I|4.6|l . are con- 
venient for numerical calculations of the form factor at 
non-zero q 2 . However, when q — > 0, the nucleon propa- 
gators in the second and third terms of the CIA result 
(|4.8|1 develop singularities that cancel, leading to terms 
involving the derivatives of the two-body kernel. The Ap- 
pendix evaluates these diagrams in the q — ► limit, and 
gives an algebraic demonstration that the singularities 



cancel. The work is somewhat lengthy, and already con- 
tained implicitly in the proof of gauge invariance. A sim- 
ilar (but algebraically different) demonstration of charge 
conservation in the two-body case was already given in 
Ref. [13 ■ 

We now turn to the expressions for the breakup cur- 
rent. 

C. The three-body breakup current 

The three-body breakup current was shown in Fig. 1261 
The final state in all the diagrams is antisymmetric. As 
discussed in Sec. Ill Bl this implies that the each diagram 
is multipied by the projection operator A 3 [defined in 
Eq. Q2.3|) ], Hence the diagrams all have the form given 
explicitly in Eq. I|2.5|) . 

The symmetrized RIA diagrams Fig. |2U (a) are 



M _ 
RIA — 



-A 3 < U 7 (p 3 , A 3 ) j^{p 3l p 3 ) ) *AiA 2 7'(Pl'P2,P 3 )+C'" 7 (P2,A 2 ) j%{p2,P 2 ) i ^x 1 \ 3 y{pi,P3,P 2 ) 



-^ 3 (l + C^ 2 3j M 7 (P3,A 3 ) j^{p 3 ,p 3 ) * Al A 27 'bl,P 2 ,P 3 ) 
\ / L J f y r y' 

-2A 3 U 7 (p 3 ,X 3 ) j^(f> 3 ,P 3 ) *A!A 27 '(Pl,P 2 ,P 3 ) • 
L J 77 ' 



(4.10) 
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Similarly, the interaction current diagrams, Fig. ( c ) and (d) are 

m d 3 fc 2 



Jf = -A3 U 7 (P3, A 3 ) 



E 2 (2tt) : 



1 + 2(V 12 



*AiA' 7 '(Pl,fe,A;3) 



(4.11) 



and the final state interaction (FSI) terms, Fig.[2j](b), arise from the core diagrams, and parallel those already given 
for the form factor in Eqs. I|4.8|l and (|4.9H 

f f m 2 d 3 k 1 d 3 k 2 x - 

FSI ~ E lE2 (2nf ^ 

H1H2 

{r Xl \ 2 \ 3 ; t _ ll ^ 2a ii{pi,P2,P3; h, k 2 ,k 3 ) G a » a >(k%) [1 + 2CPi2]j£, a (k 3 h ,k 3 ) ^W^fci, k 2 , k 3 ) 

+T\ 1 \ 2 \ 3 ^ ll 3> a (Pl,P2,P3; kl,k^,k 3 ) Gp>/3(k£)j% 1 (kZ,k 2 ) [1 + 2CPi2] U 1 {k 2 ,^ 2 ) ®viV2<x(fa,k 2 , k 3 ) 

+T\ 1 \ 2 \ 3 - P , UL2a (pi,P2,P3'i h,k 2 , k 3 ) w 7 (fc 2 ,/i 2 ) [1 + 2CP12] tfp,(k 2 , k^) Gp l p(k^)T t _ ll p a {k u fc^,*;^) j 



>; ' / / / T AiA 2 A 3 ;m^a''(Pl J P2,P3;fcl,fc 2 ,fc 3 ) G a " a '(fc 3 )[l + 2CPl2] 



X ^ka'; W ,«A fe. k 3 ) [1 + 2 CPl2] *„ lW ,a(fcl, *2, **) : 



(4.12) 



where T , A 1 A 2 A 3 :^i^ 2 a(pi,P2,P3; fci, fc 2 , fo) is the symmetrized three-body scattering amplitude with the final state on- 
shell. The off-shell unsymmetrized three-body scattering amplitude satisfies the equation shown in Fig. [7] This 
is 



1 E 

T 22\ 1 \ 2 l3;» lfJ , 2a (PuP2>P3}h,k2,k3) = g — (27r) 3 <5 3 (pi - k^Sx^M^ x 2 0,» 2 c(P2, k 2 ] P23) 



-2C 



E: 



d k 

7^)3 M 22 A 2 /3^a'(P2, fc 2 5 P2s)G a ' a " (K) 7 ^ ^ A lQ ''; MlM2 a (&2 . Pi . fc 3 ! fc l> fc 2, fo) ■ 



(4.13) 



The antisymmetrized amplitude is obtained by applying ^3 to both the initial and final state, and illustrated in 
Eq. H3J. 

I 



V. CONCLUSIONS 

Using diagrammatic techniques, we have derived a 
three-nucleon current consistent with the three-body 
spectator equations, and have shown explicitly that it 
is conserved. We obtain the current for elastic scat- 
tering, shown in Fig. |2U| and Eqs. (|4.8|) and (|4.9|l . and 
for the three-body breakup reaction, shown Fig. 1261 and 
Eqs. H4.10|) . H4.11[) . and 14.12(1 . The appendix also shows 
explicitly that this current conserves the charge of the 
bound state, even with an arbitrary choice of electro- 
magnetic form factors and in the presence of energy de- 
pendent interactions. 

Our results show that the spectator current will be free 
of any double counting if we always choose to keep the 
spectator on-shell. This simple rule leads to an organi- 
zational principal that also resolves some ambiguities in 
the choice of diagrams that might otherwise be present. 

Our derivation and discussion relies heavily on the 
beautiful method of gaugin g equa tions, developed by 
Kvinikhidze and Blankleider |13l[l7| . As already noted in 
Sec. IIII Bl the normalization of the two-body and three- 
body scattering amplitude used by KB is different from 
ours, giving different weights to various diagrams, but 



the total results arc the same. 



Another difference, also mentioned in Sec. IIII Bl is that 
KB do not rearrange their amplitudes so that the specta- 
tor is always on-shell. This rearrangement has some ad- 
vantages; it not only makes the equations more tractable, 
but it also displays how the spectator theory avoids the 
double counting problem (recall Fig. l28[) . the importance 
of which has been emphasized by Kvinikhidze and Blan- 
kleider. 

This rearrangement also leads to a different interpreta- 
tion of the diagrams in the theory. For example, Eq. (60) 
of KB includes RIA terms in which the photon couples to 
all three of the final nucleons, while our RIA terms, given 
in Fig. (fii) and (a 2 ), include only couplings to the 
particles in the final state interacting pair, with no term 
describing coupling to the final state spectator (both ap- 
proaches require the results be antisymmetrized). 

At first glance it might seem that our result cannot 
be correct - surely the photon must couple to all of the 
outgoing nucleons. But our total result does include cou- 
pling to the final state spectator as part of the FSI term 
discussed in Fig. [2HI (since this term arose from the re- 
arrangement, it is not present in KB). In Fig. [221 the 
off-shell equations of Fig. [§| are used to show that the 
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< a > (b) (c) 



FIG. 32: (Color online) Diagram (a) is the coupling to the specta- 
tor particle included by KB as part of their RIA contribution. This 
diagram is missing from the RIA part of our result shown Fig. 1261 a). 
but can be transformed using the off-shell equations shown in Fig. 
|5]into diagram (c), included as part of the FSI shown in Fig. l26f b1, 
The KB formalism does not include a term like (c) , so both results 
are identical. 



coupling to the final off-shell spectator given in KB [dia- 
gram (a)] can be transformed into diagram (b) using the 
equation of Fig.[5f c) which can be transformed using Fig. 
EJb) into the diagram of Fig. [23 ln the KB approach, 
this term is part of the RIA, but in our approach it is part 
of the FSI. This shows that care must be taken in com- 
paring the results of KB with ours, but if this comparison 
is done carefully they appear to be in a gree ment. 

Independently, Adam and Van Orden [Tj| have derived 
the same current operators using a completely algebraic 
approach along the lines developed in Ref. 01 ■ They ob- 
tain the same results, with the same weight factors. Their 
result naturally leads to an ordering in which the spec- 
tator is always on shell, in full agreement with us. These 
two approaches reinforce and complement each other. 



We plan to use the results of this paper to calculate 
the high energy breakup processes recently measured at 
JLab. 
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APPENDIX A: DEMONSTRATION OF CHARGE 
CONSERVATION 

The explicit demonstration that charge is conserved 
begins by expanding the diagramsll9l fci ) and (C2) using 
Fig. The result is shown in Fig. The 8 diagrams 
are organized into 4 pairs, with an integral over a three 
momentum k2. All the contributions from these 8 dia- 
grams can be written 



J cl 



+c2 



3e 



^X.X'a'ikl,^,^) [1 + 2CP12] 



(Al) 



where the integrand is the sum of two contributions that each become singular as q — > 

3c = 3y a'\^ a (ki, k 2 k 3 , k' 2 k' 3 ) — jc(ki, k 2 P 23 , ^2-^23) 



id 3 k 7 



VxM',Y(klk 2 + q;P^)G Y7 (Pi 3 - k 2 ) 



m- $ 2 — 



E 2 (2tt) 3 ^ 

+y\' 2 ,a'-i3\ 1 '{k' 2 \k 2 ;P 23 )G 1 - 1 {P 23 - k 2 ) A(k 2 ) j^(k 2 ,k 2 - q) 



■,f N {k 2 + q, k 2 ) A(fc 2 )J V 0tT> x 2<a (k 2 , k' 2 ; P^) 



1 



m- #2+ 



0'0 



V^ rM Ak2-q,k 2 ;P 23 )\ , (A2) 



with the single nucleon current and A(k 2 ) = (m+ j( 2 )/(2m) the positive energy projection operator for a particle 
with four-momentum k 2 . The integrals in 1A2|1 are both given in terms of the on-shell four momentum k 2 = {E 2 , k 2 }. 
The two-body kernels, V, conserve four-momentum, and are expressed in terms of independent variables, which we 
choose to be the initial and final momenta of particle 2 (usually on-shell) and the total momentum of the two-body 
system, so that generically 



V = V(k' 2 ,k' 3 ;k 2 ,k 3 ) = V(k 2 ,k 2 ;P 23 ) 



(A3) 
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+3: 




+12' 





FIG. 33: (Color online) Diagrammatic representation of expansion of diagrams 1191 I'd ) and (C2). The first two rows use Fig, 1311 and the 
second two rows use the on-shell version of the same figure again (to replace the amplitudes enclosed in the dashed boxes). The resulting 
8 diagrams can be organized into 4 pairs, which members of each pair arranged above and below each other on lines 3 and 4 of the figure 
[diagrams (ai) and (a2) are an example]. 



with P 2 3 = P - fa = fa + fa = k' 2 + fa. 

To calculate the q — > limit of l|A2jl . we follow the method developed in Ref. Since the full current is gauge 
invariant, we know that the singular terms must vanish as q — > 0, so we can find the q — ^ limit of ljA2fl by contracting 
q^ into both sides, expanding about q^ = 0, and extracting the coefficient of the linear term. Using the Ward identity 
for the nucleon current 



and, for fa on-shell, 



*) = (m-#)-(m- f) 



(to- fa)A(fa) = 0, 



Eq. HA2(1 reduces to 



TO (Pfa 

E 2 (2tt) : 



23 ) 



(k' 2 \ fa + q; P^)G^{P^ - fa) k{fa) >pVp (fa,k' 2 ;P 
-V x > a ^, Y (k'lfa;P^)G^(P^ - fa)K(fa)p ip Vp, rM , a {fa - q, k' 2 ;P^ 3 ) 



(A4) 



(A5) 



md?fa v - 
a" 



E 2 (2tt) 



<) 



dk 



2p 



■V K , a/]K< -y'(k 2 \fa;P^) 



+V\' 2 , a ';\q, J '(k 2 ,fa;P 23 )G Yl {P 23 - fa) 



G ~f'i( p 23 ~ fc 2) V X n n] x 2 , a (fa,k' 2 ;P. 
-V\'' 1 -\ 2 , a {fa,k' 2 ,P 23 ) 



23 ' 



2/i 



-Va' a':\''., 1 '{k 2 , k 2 ;P 23 ) 



d 



dP 



— Gy 7 (P 2 ' 3 - fc 2 ) 



23 p 



^" 7 ;A 2 , Q (fc2,fc 2 ;P 2 '3) +0( 9 2 ). (A6) 



Use the shorthand notation 
d 



dk 



2/' 



2 /i 



dG(fc 3 ) 
dk 3fJ , 



G^ifa) 



(A7) 
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+ 4 




(D (2) 



FIG. 34: (Color online) The second two-body scattering amplitude in Fig. (c) can be replaced by its integral equation, and further 
simplified using the on-shell version of|5](b). This form of the vertex function that is useful in the reduction of diagrams I19l fei 1 and (e2). 




= 12 



+12 



FIG. 35: (Color online) Diagrammatic representation of expansion of diagrams 1191 fei ) and (e2). The first two rows use Fig. I31I and the 
second two rows use Fig. I34I (to replace the amplitudes enclosed in the dashed boxes). The resulting 8 diagrams are labeled by both the 
parent diagram (a-d in the first two rows) and the term (1 or 2) in Fig. 1341 from which they originate. They are organized into 4 pairs, 
with members of each pair arranged above and below each other on lines 3 and 4 of the figure [diagrams (a2) and (c2) are an example]. 



where, in general, <5f refers to the derivative with respect to the initial k 2 momentum, and 8^ with respect to the final 
k 2 momentum. With this notation 



lim j% 



J {^, Q ';A^y(fc^fc 2 ;^3)G^ 7 (fc 3 )^, 7;A2 , Q (fc 2 ,fci;^3) 

5iV\> 2 , a > ; \» t y(k'2, k 2 ; P 23 ) G 7 ' 7 (fc 3 ) Vx^ 7; A2,a( fc 2, k' 2 ; P 23 ) 



(A8) 



Inserting this back into the original expression I)A1(I . and using the wave equation [ (|4.6[1 with the second particle 
on-shell], gives 



lim , = 3e 



^#^(2^)6 X! 1 r AiA 27 '( fc l. fc 2,fc3) G^ 7 (fc 3 )r AlA27 (fci,/c2,fc3) 



A1A2 



(fc 2 ,fc 2 ;P23)*A 1 A 2 a(fcl,fc2,fc 3 ) 

^2 ( Z7T ) w 

A 2 

/TYL ]\! — 
g , *^w(kl,k' 2 , k' 3 ) ^Vx' 2 ,a';A 2 ,a(fe 2l fc 2 ;P 23 ) [1 + ^Pu] * Al A 2 a (fcl , k 2 , k 3 ) 



(A9) 



Note that the relative signs of the terms change because of the sign in Eq. (|4.6|) . 

We now turn to diagrams Fig.^|(ei) and (e2). The object here is to isolate the electromagnetic coupling in a loop 
as we did above, but because the loop now involves a spectator, we need to pull out two interactions. In preparation, 
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we note that the Fig. |5| (c) can be written in an alternative form that involves the kernel, as shown in Fig. [M] The 
algebraic form of this equation is 

f m 2 d 3 k' d 3 k' sr-^ 

r Al A 2 a(fci,fe,fc 3 ) = 2CV 12 / 1 6 2 22 M W l7 '(ki, k[;P 13 ) G 7 > 7 (P 13 - k[) 

J £j 1 £j 2 \ A7T ) w w 

A 1 A 2 

xV X2 ^ Ka ,(k 2l k' 2 ;P^) [1 + 2CP 12 ] y x[ x> a '(k[, k' 2 , fc 3 ) , (A10) 

where Pi 3 = P — fc 2 and P 23 = P — k[. 

Using this substitution, the diagrams Fig. [Wl (ei ) and (e2) can be written as shown in Fig.OHl These 8 diagrams 
collect together into 



J el+e2 ~ izc 



m 4 d 3 k'{d 3 k[d 3 k!{d 3 k 2 



E»E> X E<>E< % 



j^yT- *K*M k l *£> *s) [1 + 2CP12] 



X^S'U.J^^fcs. *1*2*3) [1 + 2CP12] * Al A a a(fci, 4, 4) 

where the common internal loop is again the sum of two (canceling) singular terms 

/TTl k'2 
2jT (27r) 3 

1 



(All) 



(4> fe 2 + g;p23) 



m- # 2 - ^ 

+^ > a':.-.-l/r?./,,:P^i A(fc 2 ) j£(fc 2 ,fc 2 -c,i 
where 



jft(fc 2 +<?,fc 2 ) A(fc 2 ) 



Vfl, 7; A 2 , Q (fc2, & 2 ; ^23) CA' l7 ',A 17 (fci , fei; P13 - q) 



1 



/3'/3 



^, 7 ;Aa,a(fc2-5,fc 2 ,P^)OA' (k'{,k[;P 13 ) , (A12) 



Ai y iAl7 (fc;',fe;;Pi 3 ) = G 7 ^(P 3 - M') M AX)Alw (fcf,fe;;Pi 3 )G U)7 (Pi3 - fej) 



(A13) 



and P 23 = P — k[, P 23 = P' — k'{ , P' = P + q, and P13 = P' — fc 2 . Following the same arguments that lead from 
Eq. (|A2[I to JXSJl, we obtain 



A" 7 ;A 2 ,a(^2, fc 2 i ^23) 



- <5f y^,a';A» 7 '( fc 2! fc 2 ; P 2 " 3 ) Ai y iAl7 (fe", -P^) ^A 2 ', 7 ;A 2 ,a(fc2, fc 2 ! P 23 ) 

-^,a';A»,y(fc 2l fe2;P^)C (fci',A4;i^) ^FA», 7; A 2 ,a(fe,A; 2 ;P^)} , (A14) 

where (J^V and V have been previously defined, and 

A»0 KYMl (k';,k[;P 13 ) = G%,(P 13 - k'D M AyiAlW (^' ! fci;P 13 )G £l)7 (Pi3 - 

+ Gy w /(Pia - k'l) M Kui , Mu {k'(, k[;P 13 ) G£ 7 (P 13 - k[) 
+ G 7 ^(P 13 - fc?) S^My iU , M Jk';, k[;P 13 ) G W7 (P 13 - fci) , 

with 

d ^ 



(A15) 



8>iM{k'l,k , 1 ;Pv. 



ap 



13 



M(k'{,k[;P 13 ) 



(A16) 



the derivative with respect to the total two-body momentum. Inserting l|A14|l into IjAllfl . and using the replacement 
(IA15|) and the wave equations l|4.6|l and (|A10|) . gives an intermediate result 



f f m 2 d 3 fcxd 3 fc 2 ^ J- 



fe, * 3 ) g^ (fe 3 ) 4CP12 r AlAa7 (fci, fc 2 , fc 3 ) 



+ 



E ^Aia'^l.^. fc 3) ^.a-A^a^, k 2 :P 23 ) M> Al A 2 a(fcl , *2, fcs) | • (A17) 



18 



where 

<*T>';A 2 , Q ( fc 2, k 2 - P 23 ) - [1 + 2(V 12 ] S?Vy 2 , a , Mia (k 2 , k 2 ;P 23 ) 2CV 12 + 2CJ> 12 5»V x ^ a , M , a (k' 2 , k 2 ; P 23 ) [1 + 2(V 12 ] 

-2CV 12 SpMy 2i a'-x 2!a (k 2 , k 2 ;P 32 ) 2(P 12 . (A18) 

The evaluation of the derivative of M can be carried out using the fact that M is an infinite series of interactions: 

M = V - VGV + VGVGV , (A19) 

and recalling that k 3 = P 23 - k 2 , so that dG(k 3 )/d(P 23 ) fl = dG(k 3 )/d(k 3 ) p = G^ 

6 P M = 6 P V - {S P V)GV - VG^V - VG(S P V) + (S P V)GVGV + VG^VGV + VG(S P V)GV + ■■■ 

= 6 p V - (S p V)GM - MG(5 p V) - MCM + MG(S p V)GM . (A20) 

Therefore, the 8 P M term in ilA18(l becomes 

Ai 



E „ {2n 2 )3 E SpVy 2 a',x»p>(k' 2 , *?; ^23) G , {k'£)M x „ 0Ma {k'l fc 2 ; P 23 ) 
- J E^^f ^ M \' 2a >,\'if3>{k' 2 , k 2 ;P 23 ) Gp>p(k3)8 P Vx»/3,\ 2a {k' 2 , k 2 ;P 23 ) 



E „ (2 ^ 3 E M y 2 »',y^'(k' 2 , k'i;P 23 ) Gfo(k%)Mx, PtXaa (kZ, k 2 ;P 23 ) 

xG p , p (k' 3 ") M W n pMa (k' 2 \ k 2 ,P 23 ) . (A21) 



Using this expression and the original form of the wave equation, (|2.9[) . the term involving the derivative of M reduces 
to 



lim j£ , 



= 3e // ^ F k (OTA6 E I r AiA 27 '(^1^2,fc3) G^ 7 (fc 3 )r AlA27 (fci,/c 2 ,A:3) 



Me term J J E\E 2 (27r) 6 

Ai A2 



/md?k' 1 
E , (27r) 2 3 E *AiAia'(*l. 4, ^3) [1 + 2CP12] S P V XWMa (k 2 , k 2 -P 23 ) [1 + 2C^i 2 ] * Al A 2Q (fcl, k 2 , k 3 ) \ . (A22) 

^ A~ J 



Collecting all terms from Eqs. HA9j) . (|A17fl . and (|A22jl gives 



lim J^j_„ = 3e 



// %t 2 [tf ^ {rA 1 A 2 y(fci,fc 2 ,fe3)G^ 7 (fe3)[2 + 4CPi 2 ]r AlA27 (fc 1) fc 2 ,fe 3 ) 

A1A2 k 

f m d^h' ^ - 1 

+ y -^ 2 ^E^w( fc i< fc 2 I fc 3 )[l + ^^ (A23) 

A 2 ' 

where 

(A^ +e = (5f + ^-^)^ )a , A2 , a (fc 2 ,fe 2 ;P 23 ). (A24) 

We next look at the q — » limit of the interaction current terms, shown in Figs. EH (f) and written in Eq. (|4.9|l . 
Rewriting the WT identity for the interaction current, Eq. I|3.7|) . in terms of the independent variables gives 

q^Vf = V(k' 2 , k 2 -P 23 ) - V(k' 2 , k 2 ;P 23 + g) + V(k' 2 -,k 2 ; P 23 ) - V(k' 2 , fc+; P 23 + q) . (A25) 

Expanding l|A25|) in powers of and equating the coefficient of the term linear in g M , gives 

l ™ V i%,*UXr,J k ^ k * P ™) = -{s^+S^ + 2Sp}Vy 2ia , X2 . a (k 2 ,k 2] P 23 ). (A26) 



19 



Substituting this into i|4.9|) gives a result with the same form as the second term in (|A23|) . and combining this with 
(|A24() gives a total contribution of — 35 pV. Combining this with the contributions from diagramsll9l fb) and (d), and 
the other terms from (|A23|1 . gives the total result 



•Cal(O) 
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wr d A kid A k 2 
E X E 2 {2%y 



J2 r Al A 2 y (fci, k 2 , k 3 ) G^ 7 (*a) [1 + 20Pi 2 ] r AlAa7 (fci, k 2 , k 3 



-9e 



Ai A2 

m 3 d 3 kid 3 k2d 3 k' 2 
E 1 E 2 E 2 (2irf 



*A 1 ^o'(*l, k' 3 )[l + 2(V 12 ](6 p V)[l + 2(V 12 ] *A 1 A 2 a(fcl, fc 2 , A* 



lim 3eu(P', A' 



Ai A2 A 2 



2M B 



u(P, A) = 3e 



Mb' 



(A27) 



where the last line uses the fact that the total charge of a bound state of three identical particles of charge e is 3e 
(because isospin has been ignored) . Hence charge is conserved if the normalization of the wave function is 



1 = 3 
Noting that 



PA* f)C 

r — ^— - [1 + 2CVi 2 ] r 

M B dP^ 



P dV 

* [1 + 2^12] ^^[1 + 2^2] * 



P„ dG 
Mb dP^ 



= 2Mb 



dG 
W 1 



(A28) 



(A29) 



and similarly for V, we recover a normalization condition equal to 2Mb times that originally derived in Ref. |ll| . 
This difference is due to the fact that the all spinors in Ref. [Il| the were normalized to 2m (or 2Mb for the bound 
state), so our result agrees with Ref. pl|. completing our demonstration. 
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